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Abstract 

The Standard Model predicts that the Higgs boson couples to the fermions in the mass eigen¬ 
states. We consider the effects of lepton flavor violating (LFV) Z boson couplings in the case where 
the Higgs boson has flavor non-diagonal Yukawa interactions with the muon and the tau lepton 
generated from physics beyond the Standard Model. We list the formulae of the couplings of the 
effective interactions among the r lepton, the muon and the Z boson. Using these formulae, we 
calculate the branching fractions of various leptonic and hadronic LFV r decays, and the LFV Z 
boson decay: Z rp,. Although the Z-boson contributions to LFV tau decays cannot be ignored 
in terms of the counting of operator dimensions or chirality flipptings, it turns out that they are 
not very significant for r —>■ 3/i and t ^ pp decays. We also calculate the branching fractions of 
the processes, r —>■ pn, r —>■ prj^'l and t — >■ pai, which are dominated by the Z-boson exchanges 
due to the spin and the parity of the hadrons. 
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1 Introduction 


The Higgs boson has been discovered at LHC mm in 2012 with mass around 125 GeV mm, and all 
particles in the Standard Model (SM) are now discovered. Nevertheless, the nature of the Higgs field 
and its potential is still unknown. It is quite natural to expect that there is some larger framework, 
physics beyond the standard model (BSM), that provides us with better understanding of the nature 
of the Higgs boson. 

The lepton flavor violation (LFV) is one of the clear signals of BSM. Although the observation 
of neutrino oscillation phenomena (see Ref. [5] and references therein) implies that lepton flavor is 
not conserved in the neutrino sector, simple incorporation of the neutrino mass in the SM does not 
result in the charged lepton flavor violation at the observable level laiz!. Since the Higgs field is the 
origin of the flavor structure in the SM, it is quite conceivable that BSM hidden behind the Higgs 
mechanism induces the LFV processes with their rates much larger than the ones predicted from the 
neutrino mixings. Indeed, it has been shown that various scenarios of BSM predict large branching 
ratios of LFV processes, for example in supersymmetric models [El El uni dn na [i3], multi-Higgs 
doublet models [H], the littlest Higgs model with T-parity dSlIIlKIZj, and the SM extended by extra 
dimension [I81II9]. 

After the recent discovery of the Higgs boson, there have been many studies on its properties such 
as spin-parity and couplings Oil!. One of the non-trivial predictions of the SM is that the Higgs 
boson couples to fermions in the mass eigenstates [20]. This prediction is not necessarily true in 
BSM physics, and thus searching for flavor off-diagonal couplings opens up the opportunities to probe 
the physics behind the Higgs mechanism. The possibility of LFV couplings of the Higgs boson has 
been originally studied by Bjorken and Weinberg m (see also mm), and after the start of the 
LHC experiments the model has been paid a renewed attention mmmm- In these works LFV 
processes induced by the off-diagonal Higgs couplings have been discussed, especially fi —)• ey, r —)■ /ry 
and r —ey decays. Bjorken and Weinberg have pointed out that two-loop diagrams provide the main 
contribution to these processes rather than one-loop diagrams due to the chirality structure. It has 
later been estimated that the two-loop contributions are a factor of a few larger than the one-loop 
ones [Ml[^ . 

LFV processes have been searched for in the decays of the muon and the tau lepton. The current 
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experimental limits on the branching fractions are i3(/i —>• 67) < 5.7 x 10“^^ [28], —)• 3e) < 

1.0 X 10“^^ [21] and B{^Au —>• eAu) < 7 x 10“^^ [301 for muons. For LFV tan decays, there have 
been searches for decay modes such as r —)• /ry, r —>■ 3 leptons, and r —)• ;U + hadrons m, and the 
upper bounds have been obtained as 0(10“®). In particular, the BaBar and Belle experiments set 
a limit B{t —>• /ry) < 4.4 x 10“® [321133] . Although the upper bounds look much stronger for muon 
decays, tau decays may be more important if the origin of LFV is related to the Higgs mechanism 
as the Higgs field couples strongly to matter in the third generation. In future, the sensitivity to the 
branching ratios of the LFV tau decays such as r —)• /ry and r —?■ 3/i will be improved to at the level 
of 10“^®“^°^ |3l| at the Belle II experiment. 

In addition to the B factory experiments, there have been searches for the LFV Higgs decays at the 
LHC experiments. Recently, the CMS collaboration has reported the upper limit of the LFV Higgs 
decay into the tau and the muon using 19.7 fb“^ of ^/s = 8 TeV data set taken in 2012 [35|, and this 
provides us with the upper limit on LFV Yukawa couplings of 0(10“^) [241135] . Interestingly, the 
analysis of CMS has also reported a 2.4 a excess in the H ^ Tfi decay mode [35] . 

In the literature on the LFV tau decays through the non-diagonal Yukawa couplings, one-loop 
and two-loop photon mediated diagrams and the tree-level Higgs exchange diagrams are evaluated 
for T —)• |Uy, r —)• 3/U and r —)• /r -|- hadron decay processes [SIES]. According to their formulae the 
sensitivities to the LFV couplings at the future B factory experiments are found to be similar to those 
at the LHC experiments. 

In this paper, we discuss the Z boson mediated contributions to the LFV tau decays in the 
model with the flavor off-diagonal Higgs interactions. Such contributions have not been calculated in 
literature although, based on the dimensional analysis and the chirality structure, the contributions of 
the Z boson to effective LFV four-fermion operators are expected to be comparable to other diagrams, 
such as photon mediated diagrams. Moreover, in r decays there are hnal states which are allowed 
only through the Z boson. We first list the coefficients of effective T-fi-Z interactions for both the 
monopole and the dipole types. Although the dimension of the dipole operator is higher than that of 
the monopole type, one cannot ignore the dipole one since the contribution to the decay amplitude 
is comparable considering the chirality structure. Using these coefficients, we evaluate the LFV tau 
decays as well as LFV Z decays. It turns out that the Z boson mediated contributions to the tau 
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decays are numerically sub-dominant when there are contribution from photon mediated diagrams such 
as r —3^ and t ^ decays. For the hadronic decays into pseudoscalar and axial vector mesons, 
the Z boson contributions are dominant, and the dipole operators are found to be numerically less 
important compared to the monopole ones. 

In Sec. [21 we discuss the model of the LFV Yukawa couplings and gauge invariant higher dimen¬ 
sional operators which generate the couplings. In Sec. [3] we list the coefficients of the effective r-p-Z 
interaction obtained by evaluating the one-loop Z penguin diagrams, then by using the coefficients 
we evaluate the LFV Z decay. By using the effective Lagrangian, we calculate the branching ratios of 
LFV tau decays in Sec. jH We summarize our result in Sec. [5l 

2 Lepton Flavor Violating Higgs couplings 

Although the couplings of the SM Higgs boson are predicted to be flavor diagonal, the BSM contribu¬ 
tion represented by higher dimensional operators can modify the prediction. In this section, we review 
how such flavor off-diagonal Yukawa couplings are generated from dimension-six operators. 

One of such operators is LeR(p{(p^(p), where L and p are left-handed lepton doublets and the Higgs 
doublet, respectively, and cr is right-handed lepton fields. The Lagrangian which is responsible for 
the Yukawa interaction of lepton fields is 

^Yukawa = -VijL^e^R^ “ V) + H.C., (2.1) 

where yij and Cij are dimensionless coupling constants. Indices i,j = 1,2,3 indicate the generation 
of the lepton fields, and A is an energy scale for the normalization. The Yukawa terms Tvukawa are 
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expanded around the VEV of the Higgs boson as follows: 


■^Yukawa 


rriij 


Yi 




- iYijel&j^(j )2 + H.c. 


V2A2 


Cije\ej^ 


(f> + 2iH(f)2 + (/’i) 


+ H.C. + 








( 2 . 2 ) 

(2.3) 

(2.4) 


where in a general gauge the Higgs doublet is expanded around its vacuum expectation value (VEV), 
V, as (f = {4>'^,{v + H + i(j) 2 )IV2)'^. The charged and the neutral Nambu-Goldstone bosons are 
denoted by 4>'^ and (/>2, respectively. One can see that the factor of three in Eq. (|2.4p prevents us 
from simultaneously diagonalizing the mass and the Yukawa matrices unless Cij = 0 or Cij oc yij. In 
the mass eigenstate of lepton fields with mass matrix m = diag(me, m,-), the Yukawa matrix Y 
has, thus in general, non-vanishing off-diagonal elements. In addition, there are new flavor violating 
interactions such as the fourth term in Eq. O, which is necessary to maintain gauge invariance. 

In this work, we consider the effects of r-fi-H interaction, and we assume other off-diagonal elements 
are zero for simplicity. In the mass basis of the charged leptons, the Lagrangian is given by 


^Yuakwa = -mildR “ YiIlIrH - fj, {Y^^Pr + Y*^Pl) tH + H.c. H -, (/ = e, r), (2.5) 

where Pr = (1 -|- 75)/2 and Pl = (1 — 75 )/2. We also assume that the diagonal components of the 
Yukawa matrix is SM-like, i.e., Y) = raijv. 

These LEV Yukawa couplings are directly constrained by the searches for H ^ T\i process at the 
LHC experiments. The current experimental upper bound is B{Pl —)■ = B{H —)■ -|- 

B{H —)• r“/i+) < 1.51 X 10“^ [3S], whereas the theoretical prediction is 

B{H ^ r±/i^) = 1.2 X 10^ (|Y,^|2 + |Y^,|2) , (2.6) 

where the Higgs boson mass is taken as mu = 125 GeV and we use the SM prediction of its width 
Tj^ = 4.0 MeV [3B] at mn = 125 GeV in the evaluation. 
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3 T-^i-Z from LFV Yukawa coupling 


3.1 Calculation of vertex from one-loop amplitudes 

In the presence of the flavor violating Yukawa interactions, the effective t-^-Z couplings are induced 
by one-loop diagrams. The r —/rZ* transition amplitude is in general parametrized as follows: 

\ rrir rrir J 

{CIC{s)Pr + CE*{s)Pl) Urtliq) 

ur2iq^el{q), (3.1) 

\ rur rrir J 


where q^ is the four momentum of the off-shell Z boson, and s = q^ and = (z/ 2 )[ 7 ^, 7 ^]. The wave 
functions of the muon and the tau lepton are denoted by and Ur, respectively. The polarization 
vector of the Z boson is denoted by e^{q). We set the initial r and /r in the final state as on-shell. 

The dimensionless effective couplings ^{s), ^^{s) and Drr{s) are in general functions of s. In 

the calculation of Z boson decays and tau decays, the contribution from l{s) can be neglected since 
= 0 for an on-shell Z boson and = 0{mf) where and ruf are respectively the current 
and the mass of the light fermions in the final states. Although the terms with couplings r{s) and 
Drl{^) respectively correspond to dimension-four and five operators, the ones with Drr{s) can be 
equally important in general as we will see later. 

The one-loop level diagrams which generate the above amplitudes are shown in FIG. [H where we 
omit diagrams which are suppressed by the muon Yukawa coupling. Summing ail six diagrams we 
obtain a finite result which does not depend on a gauge parameter. 

By evaluating diagrams in FiglU we obtain 


C'f('S) — Wvi^)9v + ) 

F>l{F) = — \.Fd{^)9v + Ff){s)g^jE\ , 


(3.2) 

(3.3) 

(3.4) 

(3.5) 
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H 



FIG. 1. One-loop finite and gauge invariant set of diagrams which contribute to t-^-Z vertex. Black 
dot denotes the LFV Yukawa coupling. 


where the functions Fy^{s) and are expressed in terms of the Passarino-Veltman functions 

Bq, and Ci, {i = 1, 2,12,22,00) defined in App. as follows: 


Kis) 


— (2m^(C'i + C 2 + C 12 + C 22 ) + 4C'oo)(0, s, rn\j, m^, m^) 

+2{2ml - s, 0, ml) 

+‘im%Co{s, ml, 0, mj^, m\,ml.) + ‘im^iCo + C 2 ){s, 0, ml,m%, m|, 0) 


+ 4- 


m 


H 


2 


{Bo{ml, mlj,ml) - Bo(0, ml^,ml)) 


F^is) 


+Bo{0, mjj, ml-) - 2 Bq{s, ml, ml), 

= {2ml{-Ci -C 2 + C 12 + C 22 ) + 4C'oo)(0, s, ml, mjj, ml, ml) 

-2{mlj -h 2ml)Co{ml,s,0,mlj,ml,ml) - Am%{Co + C' 2 )(s, 0, m^, m|, 0) 

+A{-mlC 2 + mlCi 2 - m%Co + 2Coo)(s, ml, 0, mlj, ml, ml) 


(3.6) 


+ 


m 


H 


mt 


-4) {Bo{ml,mH,ml) - Bo{0,mH,ml)) 


-Bo{0,ml^,ml) - 2Bo{s,ml,ml), 


(3.7) 
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(a) (b) 


FIG. 2. One-loop amplitudes induced by both Tfj,h(j )2 and interactions in Ea. ()2.2l) which 

are denoted as black dots. 

P'nis) = -m^{Ci + 2C2 + Ci2 + C22)iO,s,ml,mH,ml,rnl) 

+2m% 0, ml, m^, 0) -b Ci{s, ml, 0, m|, ml)) , (3.8) 

Fd{s) = -ml{Ci-2C2-Ci2-C22)iO,s,ml,mH,ml,ml) 

-2ml{C2 + Ci2){s,ml,Q,mlj,m\,ml) 

+2m\ (C'i(s:0,m^,m|-,m|,0) - Ci{s,ml,^,mlj,m\,ml)) . (3.9) 

In our notation the Z-boson interaction with fermions (/ = u,d,e) which have electric charge Qf are 
given by 



- <,zh-Y''Y^PnF''-/^phz„ 

(3.10) 

9v 

= -2W{ew)Qf, 

(3.11) 

A 

_ rp^ 

(3.12) 


where 9w is the Weinberg angle and = 1/2 for up type quarks and = —1/2 for down type 
quarks and charged leptons. 

Because of the mass relation m-j- <C mz,mH, the functions F^is) and F^°'{s) can be approximated 
by expanding in terms of r,- = ml/mlj. We find Fy°'{s) and F^{s) are of 0{rl) while F^(s) is of O(rT-). 
The fact that and F^(s) are the same order stems from the specific features of this model, i.e., 

the LFV interactions are always accompanied by a chirality flipping. The leading contribution to the 
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dipole couplings comes from the diagrams (d) and (e) in Fig. [H because these diagrams are 

not suppressed by m,-. Regarding the monopole couplings every diagram is suppressed by 

nir because we need another chirality flipping in addition to the LFV couplings. Therefore, C^^{s) 
and ^(s) are at the same order in terms of r^- in the normalization defined in Eq. (13.ip . and thus 
one cannot neglect D^^{s) even though the dimension of the dipole operator is higher than that of 
the monopole one. 

The function F^{s) vanishes at the leading order of r,- due to the cancellation between the diagrams 
(d) and (e) of Fig. (T} This is again because of the specific chirality structure of this model. Once 
we hx the chirality of t~ in the external lines, the internal t~ and in diagrams (d) and (e) have 
opposite chiralities, which results in the opposite sign of the axial Z boson coupling. 

Note that there are other diagrams shown in Fig. [2] when we consider the full Lagrangian in 
Eq. ()2.2I) . Although these diagrams are not suppressed by r^, they contribute only to the couplings 
™ Eld- P'l.ip . and thus we can ignore them. 

3.2 LFV Z decay 

Using the effective couplings defined in Eq. m, we obtain the following branching fraction of Z —)• r/i. 
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1 rnz 
F^ bvr 


(|Cf (m|)p + |Cf (m|)p) 




(3.13) 




where we neglect contributions with higher order of r,-. Here we list the functions 
F^°‘{rn? 2 ) the leading order of r^-: 


Fv{rnz) 


Fv{m%) 



_ 2 (4rz 3 1) —tan ^ (\/4rz - l) + 8 (tan ^ (A/4rz - l))^ + ^Li 2 (-rz) 

_ log(j^z) ^^^2 _ 4 ,,^ _ 2 log _|_ 1) _g i) + 2 (log {rz)f 
^z 

+T “ ^ - 2rz + 2 log (rz + 1)) , 

~ ^ tan-i (V4rz - l) - 8 (^1 - (tan"! (V4rz - l))" 

+ log(j^z) + 21og (rz + 1)) + 4 -Li 2 (-rz) 

^z ’’z 

(l7rz + {Arz + 2) (log (r^))^ - 4^ 

(?'z - 2?'z + 2 log (rz + 1 )) , 

^z 

O(r-r)^_ 

— 4V4rz —1 tan"^ (V4rz — l) + — (tan“^ (-v/4rz — l))^ 

+ i5*M(log(r.) + 2) + 4, 

rz 


(3.14) 


(3.15) 

(3.16) 


(3.17) 


where Li2 is the dilogarithmic function and rz = m^lmj^. 

We list the numerical values of Fy‘^{m\) and F)(,’“(m^) as well as the effective couplings 
and D^j^{rin? 2 ) ia TAB. [H where we use the parameters listed in App. |Bj Except for Fh {m?z) which 
is of O(rT-), E’s are all comparable. However, the couplings ^(rrv^) turn out to be two orders of 
magnitude smaller than This is caused by a numerical accident in the vector coupling of 

the Z boson to the charged lepton, gy = —1/2 + 2sin^ {Ow) = —0.038. By considering the prefactor 
m^/m^ in Eq. (|3.13l) . the contributions of the dipole and the monopole interactions are comparable. 


The imaginary parts of E’s are originated from the absorptive part of the diagram in FIG. 1(c) 


This diagram contributes to Fy's at the leading order of r^, while it is not important for Fh- For Fb, 
since the real part is already 0{rT-), the imaginary part is comparable to the real one. 
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5.0 - 0.78z 


(2.3 - 0.30i) X 10-^ 


-4.8 - 0.78i 


(-2.0-h 0.35i) X 10-5 

F^{m%) 

(-8.6-M.6i) X 10-® 


-2.7 X 10-'^ 


0.84 

D^{im?^)lYr^ 

-2.7 X 10-'^ 


TABLE 1. Coefficients of gy and g\ in the effective couplings of dipole and vector operators at 
s = m\. Regarding the value of we use the functions dehned in Eq. (13.9p . 


By using these values the branching ratio is obtained as 

B{Z^ = 8.9 X 10"^°| + 7.7 x (3.18) 

The difference of the coefficients between the first and second terms is caused by parity violation in 
the weak interaction. 

4 LFV tau decays 

In this section, we discuss the LFV tau decays induced by the off-diagonal Yukawa couplings. We 
first introduce the general effective Lagrangian responsible for tau decays. We then calculate the 
coefficients of the effective operators by using the effective coupling obtained in the previous section. 
The branching fractions of the LFV tau decays are evaluated including the Z boson mediated diagrams. 
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FIG. 3. Diagrams which contribute to effective couplings of four-fermion interactions, where small 
black dot denotes the LFV Yukawa couplings and large black dot denotes loop induced LFV interac¬ 
tions. 


4.1 Effective four-fermion interactions 


In general, the effective Lagrangian of the processes r —)■ fj,X, [X = ^u/r, p, vr, ai) is given as follow^: 


£ dim .6 
eff 


/’dim .7 
^eff 


TTlr TUr 

+ ^ {9sLi^Rf^L)ifRfL)+gsRR{fLf^R)ifLfR) 

^ /=w,d,s,^ 

+9slR{RR9L){fLfR) + 9sRL{fLp,R){fRfL) 
+9vRR{RRl^9R){fRli^fR) + 9vlLiRL'y^9L){fL'y^lfL) 


+9vlLi'^Rl'"9R)ifL7t,fL) + 9vlR{rLl^PL){fR7tifR) 
+9TRRiRR^^'' 9L){fRCr^,ufL) + 9-R)ifLCr^,ufR)^ 


H.C.. 


4 E 

f=u,d,s^fi 

if) 


9R)d^{fLlufL) + CilRiTLCr^'' 9R)d^iURlufR) 


+CjRr{tr(T^^'' PL)di,{fR'^ufR) + Cijl^{TRa^’'pL)d^{fLlufL) 


+ H.C., 


(4.1) 


(4.2) 


where /’s are quarks or leptons. These effective couplings are induced from the diagrams in FIG. [3j 
Here we keep the light fermion masses at tree-level while we neglect the contributions of 0 {mf jrrir) 
in the loop diagrams. The effective couplings of the photon dipole operators come from one-loop and 

^For decays into scalar and pseudoscalar mesons, other operators such as r/rGG and r/rGG can contribute (see Celis 
et al. for detail). The operator rpGG is absent in this model since the Higgs boson do not couple to GG. 
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two-loop diagrams, and are denoted as -|- 1?^’^°°'’. They are evaluated in Refs. [TTl 

[SZlETlEaiM]. The one-loop level contributions are expressed as 


D 


D 


7, Hoop _ 

R — 

7,Hoop _ 


Qee f ^ 

2 ^tij,Yt I „ log 


2(47r)2 rrifj 

QeC ml 

2(472^^^/^- 


m: 


m 


H 


= - 1.0 X 10 “®K 


T/Z? 


TTr^r -X-log 


m: 


m 


H 


= -1.0 X io“^y; 


/XT ) 


(4.3) 

(4.4) 


where e (e > 0) is the coupling constant of the electromagnetic interactions. These couplings are 
suppressed by three chirality flippings and thus of 0(Y^) considering the normalization dehned in 
Eq. (14.ip . On the other hand, there are a class of two-loop Barr-Zee type diagrams [38] which are 
suppressed only by a tau Yukawa coupling, and they are found to be larger than the one-loop ones by 
about a factor of five, i.e., = —5.0 x YTfj,{Y*^) [2711371 [M| . 

Note here that calculated in the previous section are of 0{Yr), and another suppres¬ 

sion arises when the Z boson is replaced with its propagator. In total, the Z-boson dipole contributions 
are of 0(Y^), which is the same order as the one-loop photon dipole couplings. Therefore, we cannot 
neglect the Z-penguin diagrams considering that the dominant two-loop contribution is not much 
larger than the one-loop one. 

The effective couplings of four-fermion operators are listed below. The coefficients and 

are found to be vanishing in this model. The scalar couplings 55’s are induced by tree level diagrams 


which exchange the Higgs boson (FIG. 3(a)): 


9^SRR(L) 

9^sll{r) 


mz 


m 


YfY 
2 

H 


mz 


m 


TYfYr,, 

H 


(4.5) 

(4.6) 


where we do not neglect the Yukawa coupling of the light fermion. The vector couplings are sum of 
the contributions from the photon and the Z boson mediated diagrams as follows: 


9vhh' — d'vHH'i'y) + 9vhh'(Y), 


M) 


if) 


(4.7) 
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where H,H' = ii, L, and 


3vRR(L) (T') 

9vll{r) (T') 
dvRLi^) 
dVRRi^) 
9 vlr(^) 

SIl(Z) 


9 

oYrY; 

mjj 

47r 

9 

OiYrYr^^ 

mjj 

dvr 


otzYrY;, 

ml 

327r 


OizYrY;^ 

ml 

327r 

9 

m:; 


ml 

327r 



ml 

327r 


Wv-A)[Fv{m-FmsA\, 

(gi+s/A)[Fv{m’-FmsA\, 


{dv + 9a) [^vi^)9v 

{9i-g{)[Fm9t^ 


Fm9l] , 
Fm9l] • 


(4.8) 

(4.9) 

(4.10) 

(4.11) 

(4.12) 

(4.13) 


The couplings of dimension-seven operators are induced by dipole contribution of Z penguin diagrams 
as follows: 


= ^^^^^^Sv-ilA)yFHm + Ft{a)g‘A\, 
= ^^^^^^{4 + s’A)lFl.{0)g‘v + FU0)g‘A], 
= ^^^^^(9i+9i)iFDmgv-Fi(m]. 
= ^^^^^(gl-giHFm)gv-FmgA]- 


(4.14) 

(4.15) 

(4.16) 

(4.17) 


In the effective couplings of the photon and the Z boson, we take s = 0, which is a good approximation 
in the evaluation of tau decays. The formulae of T’s at s = 0 at the leading order of Vr are given by 


Fm = 

r2log(r^), 

l-rz 

(4.18) 

Fm = 

2 

(1 rz 2rzlog(rz)), 

l-rz 

(4.19) 

Fm = 


(4.20) 

Fm = 

6(1 rz)3 ■ 

(4.21) 
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F“(0) 

4.8 

cmiYr, 

2.2 X 10"^ 

Fm 

-4.2 


-1.9 X 10“^ 

Fm 

-8.7 X 10"® 


-2.5 X 10"'^ 

F-(0) 

0.78 

DU^)|Yr^. 

-2.5 X 10"'^ 


TABLE 2. Coefficients of gy and g\ in the effective couplings of dipole and vector operators at 
s = 0. 


Their values are listed in TAB. [2j The values are not very different from the ones for s = listed 
in TAB. [TJ 


4.2 Branching fractions of the LFV tan decays 

Below, we evaluate the various LFV tan decays which are searched at the B factory experiments. 


4.2.1 r —)• /i 

The branching ratio of r —3/U is given by 


B{t —)• /r g ) 
^ ^ 1-0 


V ’ ■ r' r' j 

^^[^(-121og (<5) -13) {\eDl\^+ \eDl\ 

(5 + 4 

T^( 9vlr~2^^slr + 9vrl~7;9^srl 
eDl 


+■ 

+■ 


|2 _i_ |^(m) |2'\'\ 
'7RR\ \'~'7LL\ j) 


12 

--5R 

6 


+ 3C!,fg + eDl(4C<S + 3CS)‘ 

.-V / ^ (ll^rk ('//,')* 1 \ f ^ (tl)* 


+ ‘^\9vRr\‘^ + ‘^\9vll\'^ + gdfl'sLLl^ 


g\\9sLL\ ^ \9 sRR\ 


(25SI + 9\7fR - + eDl (ig^yi*^ + - ^gUl^ 

'-' 7 Lr 9 vRR j 





in-7 


where is the mean life time of the tan lepton. The cut-off parameter 0 < 5 < 1 is introduced to 
avoid the singularity in the photon mediated contributions. In the numerical evaluation we assign 
6 = (2m^)2/m2. 
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This formula includes the contributions from the dimension-seven operators, which we cannot 
ignore in general especially when LFV is accompanied by chirality flipping. In the model we discuss, 
however, as we saw in the Z decays, the contribution from the dimension-seven operators is rather 
suppressed due to small gy and Fm- 

The leading contribution comes from photon dipole operators, and the contributions of four-fermion 
interactions (mainly the scalar ones) reduce the branching fraction by 9% through the interference 
terms. 


4.2.2 r —/r TT^,r],r]' 


The four fermion interactions generated by the photon or the Higgs boson exchanges or the photon 
dipole operator do not contribute to the tau decays into a pseudoscalar meson due to spin and parity. 
The leading contribution to such decay modes arises from the effective Z boson couplings. The 
branching fraction of r ^ /rvr is given by 




H 7r°) 


'Trf^ 


2567rmT 


1 - 


m: 


mi 


2 r 


(u—d) 

9vrr 


= Tt 


rnly gl{gA-9A?fl 


TTir 


2567rmT 


(u—d) 

gVRL 

,2 


+ 


(u—d) 

9vlr 


(u—d) 

9vll 


1 - 


m: 


mi 


(|C||^ + |Cf|2) 


1.9 X + 1.5 X 


(4.23) 


where gyjjffi = gyHH' ~ 9 vhh'’ (F, H' = L, R), and y = 130 MeV is the pion decay constant. For g 
and g’, the only difference from r —>• gir is the hadron matrix elements. One can obtain the amplitudes 
of T ^ gg and r —)• gg' by the replacement of 

+ iP-v.v'). ( 4 . 24 ) 


where the decay constants fp and fp are defined following the Feldmann-Kroll-Stech (FKS) mixing 
scheme ^9]. We can obtain the branching ratios oi t ^ gg and r —?■ gg' as follows: 


B{t- g-p) 


( ^4/p 
1 


2 



2 

B{t~ g~'K^), 


{P = g,g'). (4.25) 
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From the above relation, we obtain the branching ratios as follows: 


p) 

= 0.30R(r — )• p vr^) 



= 5.8 X 10"^^| Vp -h 4.5 X 

(4.26) 

B{t- p"p') 

= 0.27R(r" ^ p"7r°) 



= 5.2 X 10"^i| Vp -h 4.0 X 10"i^|y^^p. 

(4.27) 


4.2.3 T- ^ /x-pO 

The spin-parity of the rho meson, = 1“, allows the photon-exchange diagram induced from the 
dipole operator, and the dimension-seven operator and the vector operator induced from Z-penguin 
diagrams as well to contribute to the r —)• pp decay. The branching fraction is given by 


B{t- - 
rrfl 




327rm, 




2 

■l + 2p' 

1 


8 



pH')' 

| ^TRR 

P 


^(u-d) ^(u-d) 

^7r?D “r 


■'7RL 


[(4 + 2p)p] 




/ {u—d) {u—d)\ 

[PVLL +9vLR ) 


eD 


,7* 


+ 


^i{u d)* , ^iu d)* 
~r '-' 7 RL 


+ 


{L ^ i?)| 


= 5.5 X 10 


-7 


5.8 X 10 ^ITt/rP, 


(4.28) 


where p = — ^^hh' ~ ^thh'i ~ fp = 209 MeV is decay constant 

of the rho meson. Compared to Refs. [241125] . we include the contributions from four-fermion and 
dimension-seven operators. The leading contribution comes from the dipole operator of the photon 
and other contributions increase the branching ratio by 5%. 
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4.2.4 T- //-ai(1260) 


As in the case of pseudoscalar modes, axial vector mode is possible only through the Z boson. The 
branching ratio is given by 


B{t /i ai(1260)) 


rrf? 


ai 


327rm-7 




l + 2a 


8 


+ 


^{u-d) _ ^{u-d) 
^7RR ^7RL 


[(4 + 2a)a] 


+3? 


/ (u—d) (u—d) 

\9vll ~ 9vlr 


a 


{u—d)^ 

7RR 


-a 


(u—d)* 

7RL 


[3a] + (L ii) I 


= 3.5 X + 2.5 X 


(4.29) 


where a = and the decay constant of ai, /ai(= 230 MeV), is determined by assuming 

I3{t~ —)• RrCLl) = B{t~ —)• l'r‘^7r~7r~^ + Ur2TT^7r~) = 18.3% [ini ITT] . 


5 Summary 

After the discovery of the Higgs boson with mass around 125 GeV, it becomes important to check 
whether the Higgs boson has properties predicted in the SM. One important check is to see whether 
the Higgs boson couples to the mass eigenstates of fermions. For example, if there is H-r-fj, coupling 
in addition to the SM interactions, non-standard decay of the Higgs boson, H —r^u, and Z boson, 
Z —>■ r/i, as well as the various LFV tau decays can happen. 

In addition to the photon and the Higgs mediated LFV tau decays studied in Refs. [211 [IS], we 
complete the analysis by including Z boson mediated contributions. We calculate one-loop diagrams 
to generate the effective t-^-Z interaction and derive formulae as functions of momentum transfer. 
We find that at the one-loop level the results are finite and gauge invariant, even though the model 
corresponds to the addition of a higher dimensional operator to the SM. 

In terms of the counting of the V- insertions, the effective dimension-six and seven four fermion 
couplings induced from Z penguin diagrams are the same order as one-loop photon penguin diagrams 
attached to the fermion line. The contribution of Z penguin diagrams are, however, found to be 
small, because axial coupling is of 0{Y!^) and the coefficient of the vector type interaction gy is 
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numerically small in the SM. The effects of the Z boson couplings are included in the r —)■ 3// process, 
and also we derive the new formulae of the LFV tau decays into pseudoscalar and axial vector mesons 
in this model. 
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Appendix A Passarino-Veltman functions 

We define the one-loop functions Bq and Cj |l2l [33] : 

Boiplmlml) = (27r/i)2^y (AT) 

(p?,(pi= {2TTpf^ J ^ ^ (A.2) 

where = 4 — 2e, and 

Ni = — 7101 + (A.3) 

N 2 = {k + pi)'^ - ml + ie, (A.4) 

N 3 = {k + p 2 )‘^ - ml + ie. (A.5) 

The tensor integrals can be decomposed by their Lorentz structures as below, 

Cf, = Pi,^j,Ci + p 2 ,^,C 2 , (A.6) 

C/MU = g^luCoo + Pl,^,Pl,uCu + P2 ,^iP2,uC22 + (Pl,MP2,i. + P2,^iPl,y)Cl2. (A.7) 
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The explicit expression of these tensor functions are 


- 60 ( 0 , 61 , 62 ) = 
< 70 ( 01 , 02 , 03 , 61 , 62 , 63 ) = 

( 71 ( 01 , 02 , 03 , 61 , 62 , 63 ) = 

( 72 ( 01 , 02 , 03 , 61 , 62 , 63 ) = 

(700(01,02,03,61,62,63) = 

(711(01,02,03,61,62,63) = 

(722(01,02,03,61,62,63) = 

(712(01,02,03,61,62,63) = 

where A = l/e — 7 + log (Air ), 7 is the 


[ dx[-log(AB(x))+A + 21og/r], 

Jo 


1 rl-x 


[ dx [ 

Jo Jo 


dy- 


-1 


/o Jo Ac(x,y)’ 

r-l rl-x 


[ dx f 

Jo Jo 


dy 


y 


J^c{x,y)' 


p 1 pl—x 


X -y 


J^c{x,y) ’ 


rl pl—x 

dx dy 

Jo Jo 

pi pl—x 

— dx 
Jo Jo 


- ^log(Ac(x,?/)) + ^A + log/i 


dy- 


y 


10 Jo Ac(x,?/)’ 

f\ /■'-% (l-x-y)2 

yii-x-y) 


[ dx [ 

Jo Jo 


dy- 


J^c{x,y) 


Euler constant and 


(A. 8 ) 

(A.9) 

(A.IO) 

(A.ll) 

(A.12) 

(A.13) 

(A.14) 

(A.15) 


Ab{x) = (61 - (1 - x)a)x + (1 - x) 62 , (A.16) 

Ac{x,y) = -oixy - ( 022 /+ 03 x)(l - X - y) + 61 X + 62 ?/+ 63(1 - X - ?/). (A.17) 

Appendix B Input parameters 

The input parameters used in numerical evaluations are listed in TAB. [3l 
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a 

128-1 

Tz 

2.49 GeV 

9z 

0.741 

fn 

0.130 GeV 

sin^ Qw 

0.231 

fv 

-0.110 GeV 

V 

246 GeV 

f^' 

0.135 GeV 

mu 

125 GeV 

fq 

Jv 

0.108 GeV 

Tr 

4.41 X loll GeV-i 

4 

0.088 GeV 


1.78 GeV 

fp 

0.209 GeV 

rrip 

0.106 GeV 

fai 

0.230 GeV 

mz 

91.2 GeV 

nip 

0.548 GeV 

nip 

0.770 GeV 

nip! 

0.958 GeV 

rriai 

1.23 GeV 

rrin 

0.140 GeV 


TABLE 3. Input parameters. 


The decay constants of isospin-triplet hadrons are defined as 

— = (O I 17 r(p)) , (B-1) 

-iV2maJaie>^i{p) = (O | | ai(p)) , (B.2) 

-iV2mpfpe>^{p) = (O | {u-f>^u - dy^d) | p{p)) , (B.3) 

where p^ is a four momentum of hadrons. The decay constants of p and p', are defined as 

- iV2fl,p^ = (0 I + dy^y^d) | P{p)) , (B.4) 

-ifpP^ = (0 I sy^y^s I P(p)) , (B.5) 


where P = p, p'. The values in TAB. [3] are given in Ref. |44j . 
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